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ABSTRACT 

OAO 1657-415 is an eclipsing X-ray binary wind-fed pulsar that has exhibited smooth spin- 
up/spin-down episodes and has undergone several torque reversals throughout its long history 
of observation. We present a frequency history spanning nearly 19 years of observations from 
the Burst and Transient Source Experiment (CGRO/BATSE) and from the Gamma- Ray Burst 
Monitor (Fermi / GBM) . Our analysis suggests two modes of accretion: one resulting in steady 
spin-up correlated with flux during which we believe a stable accretion disk is present and one in 
which the neutron star is spinning down at a lesser rate which is uncorrelated with flux. Orbital 
elements of the pulsar system are determined at several intervals throughout this history. With 
these ephemerides, statistically significant orbital decay with a P or b = (—9.74 ± 0.78) x 10~ 8 is 
established. 

Subject headings: OAO 1657-415, Pulsar, Accreting Pulsars 



1. Introduction 

OAO 1657-415 is an accreting X-ray pulsar with 
a high mass companion. It was fir st detected by 



the C opernicus satellite in 1978 (jPolidan et al 
19781 ). The companion is believed to be 



Ofpe/WN9 supergiant (|Mason et al.ll2009l) . OAO 
1657-415 is one of 10 known eclipsing X-ray bi- 
nary pulsars. Its eclipse lasts 16% of the orbit or 
1.7 days. OAO 1657-415's moderate spin period 
(37 s) and its moderate o rbital period (10.4 days) 
([Chakrabartv et al.lll993l ) gives it a unique place- 
ment on the Corbet diagram which categorizes ac- 



creting pulsars by spin period and orbital period. 
Due to the relatively short or bital period and spin 
period of this wind-fed system iMason et al. I (l2009h 
speculated that this system may be in transition 
from wind-fed to disk-mediated accretion. The 
first part of this work establishes statistically sig- 
nificant orbital decay of the OAO 1657-415 system 
using data from the Burst and Transient Source 
Experiment (BATSE) previously on board the 
Compton Gamma- Ray Observatory (CGRO) and 
the Gamma-Ray Burst Monitor (GBM) on board 
Fermi. 
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Historically, the supergiant binaries have been 
subdivided into two groups according to the dom- 
inant mode of mass transfer: Roche lobe overflow 



or ca pture from the stellar wind ([Bildsten et al 
1997). This terminology is mismatched in that 



Roche lobe overflow describes the outflow from the 
companion while capture from the stellar wind de- 
scribes the accretion onto the compact object. An 
improved classification scheme would be to iden- 
tify two classes of high mass X-ray pulsar bina- 
ries (HMXBs): A) Roche lobe filling supergiants 
and B) underfilled Roche lobe supergiants. We 
present, in the second part of this work, evidence 
that, for the underfilled Roche lobe systems, ac- 
cretion can occur either directly from the stellar 
wind (direct wind accretion) or from an accretion 
disk created from the stellar wind (disk wind ac- 
cretion). The underfilled Roche lobe supergiants 
may then be further divided into two subclasses: 
Direct wind accretion and disk wind accretion. 

2. Background 

The different classes of HMXBs can be plotted 



in a P sp in - Porb diagram (Corbet 19861 ). or a Cor- 
bet diagram, where P sp i n is the neutron star spin 
period and P or b is the orbital period of the binary 
system (see Figure [J). 

The filled circles are Roche lobe filling high 
mass systems that have relatively short orbital pe- 
riods and short spin periods. The wind fed sys- 
tems (hollow circles) have longer orbital periods, 
avoiding Roche lobe filling, and longer spin peri- 
ods due to weaker accretion torques. The Be sys- 
tems (filled squares) display correlation between 
orbital period and spin period, possibly because 
the neutron star is further away from the donor 
star in systems with a long orbital period and it 
is sampli ng a weaker portion of the cir cumstel- 
lar disk ( W aters and van Kerkwiik 1989). Other 
systems show anomalous behavior, such a s SAX 
J2 103. 5 +4545 (jCamero Arranz et alj[2007h which 
is a Be/X-ray binary with a orbital period of 12.68 
days and a pulse period of 358 s. This unusual 
orbital and pulse period for a Be/X-ray binary 
places it among the wind fed supergiants in the 
Corbet diagram. OAO 1657-415 also appears to 
defy these groupings. With an orbital period of 
10.4 days and a pulse period of 37.1 s, it occupies 
an intermediate region between sources with mass 
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Fig. 1. — All the Galactic high mass X-ray bi- 
nary pulsars with known orbital periods plotted 
on a Corbet diagram. The hollow circles are the 
wind-fed supergiants while the filled circles are 
the Roche-lobe overflow supergiants. The filled 
squares are the transient Be/X-ray binaries while 
the hollow squares are a likely Be X-ray binary. 
The filled diamonds are the persistent Be/X-ray 
star. The filled stars are the Supergiant Fast X- 
ray Transient (SFXT). 



transfer via a stellar wind and Roche lobe overflow 
(jChakrabartv et al.lll993h . 

Identification of OAO 1657-415's companion 
proved elusive for 1 5 years after its discovery 
( Kamata et al.l Il990f) . In 1993, monitoring by 
CGRO/BATSE allowed sufficient data for the 
first orbital ephemeris to be c alculated yielding a 



mass function for the system (jChakrabartv et al 
19931 ). Assuming the pulsar has a mass of 1.4 



Mf?) and using the kno wledge of the eclipse, 



Chakrabartv et al.l (|1993l ) calculated the mass of 
the companion to be 14 M Q < M c < 18 M Q with a 
radius between 25 and 32 R Q . This yielded a blue 
supergiant of stellar class B0-B6. Attempts to op- 
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tically identify the companion failed due to a large 
column density but in 2002 a relatively bright star 
(2MASSJ17004888-4139214) was identified in the 
near-infrared that was coincident with a Chan- 
dra X-Ray Observatory p osition of OAO 1657-415 
(jChakrabartv et al.ll2002l ) and was consistent with 
its companion's spectral class. Due to the degen- 
eracy of infrared color-color diagrams, infrared 
photometry did not provide decisive results. More 
de tailed infrared spec troscopy performed in 2009 
bv lMason etaH (|2009h concluded that OAO 1657- 
415' s companion is of a spectral type Ofpe/WN9, 
which is a transitional stage between a type OB 
and a Wolf-Rayet star characterized by high mass 
loss, low wind velocities, and exposed CNO-cycle 
products. 

3. Observations 

BATSE was an all-sky monitor designed to 
study gamma-ray bursts and transient source out- 
bursts in the hard X-ray and soft gamma-ray 
bands ( Fishman and et al. 19891 ). The observa- 



tions reported here use data from the eight Large 
Area Detectors (LADs), which were 2025 cm 2 in 
area by 1.24 cm thick Na I(T1) scintillators oper- 
ated in the 20 keV-1.8 MeV band. These were lo- 
cated at the eight corners of the Compton Gamma 
Ray Observatory (CGRO). The LADs were uncol- 
limated, with each detector viewing half of the sky. 
The pulse timing analyses uses the DISCLA chan- 
nel 1 data, which consist of discriminator rates for 
the 20-50 keV band, continuously read out from 
all eight detectors with a resolution of 1.024 s. 
BATSE data used for the following analysis spans 
from MJD 48362 to 50683 (1991 April 16 - 1997 
August 23). 

GBM is an all sky monitor whose primary ob- 
jective is to extend the energy range over which 
gamma-ray bursts are obser ved in the Large Are a 
Telescope (LAT) on Fermi (|Meegan et al.ll2009h . 
GBM consists of 12 Nal detectors with a diam- 
eter of 12.7 cm and a thickness of 1.27 cm and 
two BGO detectors with a diameter and thickness 
of 12.7 cm. The Nal detectors have an energy 
range from 8 keV to 1 MeV while the BGO's ex- 
tend the energy range to 40 MeV. The current 
pulse timing analysis uses the first 3 channels of 
the CTIME data (8-12 keV, 12-25 keV, 25-50 keV) 
from the Nal detectors with a time resolution of 



0.256 s. The daily light curves are visually in- 
spected and times where the spacecraft is perform- 
ing a rapid maneuver or passing over the South 
Atlantic Anomaly are removed. Other contami- 
nating features that prevent a good background 
fit, such as gamma-ray bursts or solar flares, are 
also removed. The data from channel 0-2 (corre- 
sponding to energies between 8 and 50 keV) are 
fit with an empirical background model which is 
subtracted from the data. GBM data used for the 
following analysis spans from MJD 54690 to 55824 
(2008 August 12 - 2011 September 20). 

Swift /BAT is a hard X-ray monitor that has 
a field-of-view of 1.4 steradians and its array of 
CdZnTe detectors are sensitive in 15 - 150 keV 
range. It is a coded aperture instrument with a 
detector area of 5200 cm 2 . We use the Swift /BAT 
transient monitor results^ (15-50 keV), provided 
by the Swift /BAT team, in order to monitor the 
total flux for OAO 1657-415 from MJD 54690 to 
55824 (2008 August 12 - 2011 September 20). 

4. Data Analysis 

Because pulsars emit a large percentage of their 
radiation in periodic pulses, Fourier components 
of the source signal can be identified and ex- 
tracted from the background. For details on 



this procedure see Camero-Arranz et al. ( 2010l ) 



Finger et al 



(|l999l ). In the case of BATSE data, 
short intervals (roughly 500 s) were fit to a Fourier 
expansion in pulse phase plus a quadratic spline 
model that accounts for the combined background 
and the average source flux. In the case of GBM 
data, short intervals (again roughly 500 s) of the 
residuals from the background subtraction were 
fit to a Fourier expansion in pulse phase. OAO 
1657-415 pulse profiles, from BATSE and GBM 
data, were represented by six harmonics (n = 6) 
in this work. Additional harmonics did not result 
in statistically significant improvements in pulse 
structure. The Fourier expansion in pulse phase 
resulted in an estimated pulse profile (residual pro- 
file) for each data segment. For the purpose of cre- 
ating a preliminary phase model, the occulted in- 
tervals were removed and the remainder of the or- 
bit was initially divided into threeeequally spaced 
intervals. 



'http: //heas arc .gsfc .nasa. gov/docs/swift /results/transients/ 
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The times are barycen tered using the JPL Plan- 
etary ephemeris DE200 (IStandishlll99fih . Further 
time corrections are performed to remove the time 
delays from the orbital motion of the pulsar using 
a preliminary orbital model. This results in times 
that correspond to the pulsar emission time t em 
that are computed from barycentered times (t bc ) 
by t em = t bc — z. The line of sight delay, z, asso- 
cia ted with the binary orbit of the pulsar is given 
bv lDeeter et~aH (|l98ll ) 



z(t em ) — a x sin i [sin w (cos -E - e) + (1 
x cos oj sin E] 



(1) 



where a x is the projected semi-major axis of the 
pulsar's orbit, i is the orbit's inclination relative 
to the plane of the sky, ui is the periastron an- 
gle, and e is the orbital eccentricity. The eccentric 
anomaly, E, is related to time through Kepler's 
equation, 



E — e sin E 



■*! orb 



Tp) 



(2) 



where P or b is the orbital period and r p is the pe- 
riastron epoch, the time at which the pulsar most 
closely approaches it s companion. Initially , the or- 
bital ephemeris from iBildsten et al.1 (|1997h is used 
to correct the times for orbital motion. 

A preliminary phase model was estimated from 
historical trends in the pulse frequency. The pulse 
phase model is a monotonically increasing func- 
tion counting the number of pulses from the pul- 
sar since a defined time r referred to as the epoch. 
A simple phase model is: 



4>{tk) — 4>o + v (t k —t) + 



v o z (3) 



where v is the pulse frequency at time r and 
<p is the phase at time r. The pulse frequency, 
v(tk) is the first derivative of the phase model. In 
the simple case, v{tk) = v a , where v a is a con- 
stant and v Q = 0. As more knowledge was gained 
about the pulsar's behavior, the phase model was 
refined by adding additional terms. The long term 
phase model was determined by fitting quadratic 
splines to the corrected phases. Each unoccultcd 
interval was searched for a pulsed signal, in a nar- 
row frequ ency band, u s ing th e Y n statistic de- 
scribed in Finger et al.l (|l999l) . The search was 
performed on the first two harmonics and using 



channels 1+2 (12-50 keV) simultaneously for GBM 
and channel 1 (20-50 keV) for BATSE data. The 
preliminary frequency history of OAO 1657-415 is 
shown in Figure [2] with a long term spin- up trend 
(z> ~ 7.15 x 10~ 8 Hz day -1 ) being evident. Steady 
spin-up is evidence for a prograde accretion disk. 




Fig. 2. — Frequency history of OAO 1657-415 ex- 
tract ed from BATSE data (le ft), INTEGRAL data 
from [Barnstedt et al.l (|2008l) (center), and GBM 
data (right). 

In order to determine the most accurate phases, 
the occulted intervals were removed and the re- 
mainder of the orbits were divided into six inter- 
vals. The residual profiles were folded for each 
of th ese intervals us i ng the most recent ephemeris 
from lBildsten et ail ( 19971 ) and the best frequency 
from the above search. Although no frequency 
search was performed at this stage, the I2 statistic 
was calculated for the previously determined fre- 
quency for GBM data. Templates for the pulse 
profile were calculated by averaging the folded 
residual profiles over three orbit intervals. For 
GBM data, the 12-25 keV energy band was used. 
Long term average pulse profiles from GBM data 
in three energy bands (8-12 keV, 12-25 keV, and 
25-50 keV) are shown in Figure [3l The morphol- 
ogy of the pulse changes from a relatively smooth 
single pulse at low energy to a double peaked 
asymmetric pulse with a fast rise and slow decay 
at higher energies. There was very little evolu- 
tion of the pulse profile with luminosity or time. 
The 12-25 keV pulse profile is representative of the 
template profiles used. The template profiles were 
normalized and fit to the folded profiles along with 
a phase model. A scaling factor and phase off- 
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set was calculated from the fit. For BATSE data, 
phase offsets that could unambiguously determine 
the phase to 3<r were selected (i.e. phase offsets 
with error < 0.15). For GBM data, a cut in Y 2 
(> 8)statistic was used to select good intervals. 
The phase model was then corrected for the fitted 
phase offsets to produce the final corrected phases. 
The corrected phases and barycentered times were 




Fig. 3. — Average pulse profiles extracted from 
GBM data for lOorbits of OAO 1657-415 for 8-12 
keV (top panel), 12-25 keV (middle panel), and 
25-50 keV (bottom panel). Two cycles are shown 
for clarity. The points with error bars are approx- 
imately statistically independent 

fit to an orbital model plus a polynomial (see Fig- 
ure S]). The first term of the polynomial is the lin- 
ear trend in time expected from the phases while 
the second term represents the frequency change 
of the source over the interval. Higher order terms 
account for torque noise inherent in these systems. 



5. Timing Results 

Intervals where pulsations were clearly detected 
and where the phases did not deviate radically 
from the phase model (minimal torque noise) were 
chosen from BATSE and GBM data. Orbital fits 
to these data were performed for three orbit inter- 
vals (the same intervals in which the pulse profile 
templates were developed) from the selected data. 
Table [T] shows the intervals and results of each fit. 
A new orbital ephemeris was determined by av- 
eraging the orbital elements (with the exception 
of the orbital period, orbital period derivative and 




55780 55785 
Time (MJD) 



Fig. 4. — Black dashed curve is the pulse 
phase model minus a linear trend in emis- 
sion time for one three orbit fit. The black 
triangles are the measured phases using 
times corrected for the estimated orbit with 
the same linear trend in barycentered time 
removed. The red diamonds are the same 
phases with the same linear trend removed 
using times not orbitally corrected. The red 
curve is the best fit model of the phases 
with the same linear trend removed. The 
lower panel displays the residual to the fit 
of the orbital model. 

epoch) from each fit (See Table [2). 

Table 2: Average Orbital Elements for OAO 1657-415. 



Epoch 


T 


52298.01(79) MJD 


Orbital Period 


Port 


10.44729(21) days 


Orbital Decay 


Porb 


(-9.74 ±0.78) x UT S 


Proj. Semimajor axis 


a x sin i 


106.157± 0.083 lt-sec 


Long, of Periastron 




92.69°± 0.67° 


Eccentricity 


e 


0.1075± 0.0012 


Pulsar Mass Function 


MM) 


11.473± 0.027 M 



All the orbital epochs are plotted against time 
with a linear trend removed (see Figure [SJ. A 
quadratic fit in orbit number (see Figure [5]) 
from all the determined epochs reveals the long 
term orbital trend in the orbit of OAO 1657-415. 



The IChakrabartv et al.l (J1993J) epoch is excluded 
in this calculation because the data they used 
overlaps the data used here. Nevertheless, the 
Chakrabartv et al.l ([19931 ) epoch is consistent with 
the these results. A second order expansion of the 
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Fig. 5. — Diamonds are orbital epochs for each 
ephemeris with a linear trend removed. The line 
is a quadratic fit to the epochs. The epochs to 
the left are calculated from BATSE data while 
the epochs on the right are cal culated from GBM 
data. The epoch calculated bv lChakrabartv et al 



(|19931 ) is labeled as a red star and is excluded in 
the quadratic fit. 

epoch in orbit number, 

TV/2 = To + nP orb + -n 2 P orb P or b, (4) 

along with the quadratic fit above determines the 
orbital decay, orbital period, and epoch for this 
system (See Table [2J. Equating the quadratic 
term in the expansion of the epoch with the 
quadratic term in the fit to the epoch (Figure [5]) 
gives a P orb = (-9.74 ± 0.78) x 1(T 8 . The in- 
clusion o f the epoch calculated from INTEGRAL 
data by iBarnstedt et al. ( 20081 ) was rejected on 
the basis that it significantly reduces the quality 
of the quadratic fit (A X 2 > 20,0 00) in Figure HQ 
This suggests that the error in the Barnstcdt et al. 
(2008) epoch is significantly underestimated. 



6. Torque-Flux Correlation Results 

We compare our measured frequency rate (v) 
with Swift /B AT fluxes (Figure [6]). In order 
to determine a calibration for the BAT rates, 
observations (~ 36 ks) of OAO 1657-415 by 
RXTE were obtained on 2011 August 11. The 
RXTE/VCm (5-50 keV) spectra were fitted 



in XSPEC 12.70 with a model that includes a 
low-energy ab sorption , a p ower law, the local 
model fdcut (jTanaka 198a), an iron line and 
and a few residual lines, in XSPEC notation 
phabsx (^gaussian+power law) x fdcut. The en- 
ergy flux from 5-50 keV during these observations 
is 4.15 x 10~ 9 erg cm~ 2 s . Comparing this en- 
ergy flux and the rates from the Swift /BAT tran- 
sient monitor (15-50 keV) for approximately the 
same times, a conversion factor of 9.5133 x 10~ 8 
ergs counts -1 was determined. Rates from ob- 
servations of OAO 1657-415 by the Swift /BAT 
transient monitor were folded on OAO 1657-415's 
orbital period and converted to fluxes using the 
above conversion factor. The maximum flux dur- 
ing an orbital cycle was chosen as a proxy for 
the intrinsic flux of the source. This is because 
the orbital flux profiles suggests OAO 1657-415 
is partially obscured during much of its orbit. 
This is consistent with locally l arge column den- 



sities inferred for the source (jNaik et all [2009; 



Audlev et al.l 120061 ) . The maximum of the folded 



fluxes during the intervals in which the deriva- 
tive of the frequency was determined is plotted 
against the frequency derivative and is shown in 
Figure [5] The intervals where the frequencies were 
not monotonically increasing or decreasing were 
rejected to insure accurate determination of the 
frequency derivative. The intervals were at least 
12 days to ensure the maximum flux during an or- 
bital period could be determined. Figure [B] shows 
a strong flux-;/ correlation when v is greater than 
3 x 10~ 12 Hz s" 1 . 

7. Orbital Decay Discussion 

Orbital decay provides evidence that the sys- 
tem is evolving toward Roche lobe over-flow. Us- 
ing values for the Roche lobe radius (Rl) and 
the masses of the neutr on star and companion 
~ (|2012l ). we predict that OAO 



from Mason et al 



1657-415 companion will fill its Roche lobe in (8- 
17)xl0 4 years. This is contingent on the com- 
panion maintaining its current evolutionary pat- 
tern for that length of time. If OAO 1657-415's 
companion is experiencing the onset of a Wolf- 
Rayet stage then it is probable that it will main- 
tain its current mass loss rate and OAO 1657-415 



2 Goddard Space Flight Center. Arnaud K.A, Dorman B, 
Gordon C, HEASARC Software Development 
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v [10" Hz s~'] 



Fig. 6.— Swift /BAT flux vs. the spin-up of OAO 
1657-415 measured using GBM data. The curve 
is a model (v oc F?) that is fit to the data where 
spin-up is greater than 3 x 1CP 12 Hz s _1 . 

will evolve into a Roche lobe filling system with a 
Wolf-Rayet star given that the Wolf-Rayet phase 
lasts (3-5) xlO 5 years (jPrantzos et al.lll986l ). 

A strong stellar wind will transfer angular mo- 
mentum from the system to a halo surrounding 
the binar y system which will lead to orbital decay. 
Following IShore et al. ( 1992 ) , the system may be 
viewed as a massive star loosing material, via stel- 
lar wind, to a ring beyond the point of the sys- 
tem. If a is the semimajor axis of the orbit then 
a e is the radius of the ring of escaping material. 
In this case a £ must be larger than the L2 point 
which is roughly located at 1.2a thus a € > 1.2a. 
The ratio of escaping angular momentum per unit 
mass to total angular momentum per unit mass is 
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(M c 



M p f 



M r M n 



a(l-e 2 ) 



(5) 



where M c is the mass of the companion, M p is the 
mass of the neutron star, e is the eccentricity of 
the orbit. The predicted orbital decay from mass 
loss through a stellar wind is 



Porb ,-. „\M C +M P 

— = -(l + 3 7 ) : 

nrh 



' M c +M v 



( M a , M P 
\M C M p 



3ee 



(6) 



Assuming the eccentricity changes slowly and 
\M C \ > \M P \ and \M C /M C \ > \M P /M P \ then 



P 



orb 



P, 



orb 



-(1 + 37) 



M r 



M c + Mn 



,Mc 

'Mr 



(7) 



Using M p < 1.8 M©, and assuming M c < 18 M , 
and \M C \ > 1.0 x 10" 6 Mgyr" 1 from lMason et al 
(|2012h then 



P 



orb 



> 2.0 x 10 -6 



yr 



(8) 



where P or b is the period of the orbit. This results 
in -P or b > 5.9 x 1CP 8 . This calculation places 
a lower limit on the angular momentum trans- 
ferred from the companion to the halo via a stellar 
wi nd. Althoug h spin-orbit quadrupole coupling 



([Lai et al.l 119951) almost assuredly plays a role in 



orbital decay, the measured orbital decay for OAO 
1657-415 is consistent with the predicted orbital 
decay from stellar wind mass loss. 

8. Flux- Torque Correlation Discussion 

Flux-torque correlations are expected when an 
accretion disk is present. To test for the existence 
of an accretion disk ass uming a high M, a model o f 
disk accretion v oc F§ ( Rappaport and Joss! 1977 ) 
is fit to the data where v > 3 x 1CP 12 Hz s _1 
(see Figure [6]). Using this model, the resultant fit, 
and a range of values for the other parameters, a 
dista nce to the OAO 16 5 7-415 can be calculated. 
From Ghosh and Lamb ( 19791 ). the neutron star 
spin frequency derivative, v, may be calculated in 
terms of the moment of inertia, /, mass of the 
neutron star, M p , mass accretion rate, M p , the 
radius of the inner point of the Keplerian disk, r e , 
and a function, n(u) s ), which is of order unity and 
only depends on the fastness parameter, ui s =o'- 



v = n( Wo )(2^/)- 1 (GM p r ) 1 / 2 M p , 



(9) 



were G is t he gravitational constant. Following 
Wand ( 19961 ) the inner point of the Keplerian disk 
is given by: 



k(2GM p y^ 7 ^ 7 M p ^ 7 , 



(10) 



where y, is the magnetic dipole moment and k is 
such that r — kr a where r a is the Alfven radius 
for spherical accretion. Assuming at least partial 
threading of the accretion disk by the magneto- 
sphere, k is less than unity. 

Assuming a dipole magnetic field and that the 
gravitational energy of the accreted material is 
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converted to X-rays i.e. F — a G ^ p then, 



i> = 2.46 x 10" 



10 45 g cm 2 



1 ( R P \ 

\ 10 b cm J 



-3/7 



12/7 



( B ) 2/T f_d\W ( F \ 6/ } m 

\ 10 12 gauss J V cm / \ orgs cm 2 s ^ / 



where a is the bolometric 
the neutron star radius, d is 
the neutron star in cm, and 
nctic field of the neutron star. 



fraction, R p is 
the distance to 
B is the mag- 
By choosing a 



modest range of possible neutron star masses, 
1 < M p < 1.8, models of the equation of state 
of neutron stars determine M p , R p , and /. From 
Figure this range of neutron star masses re- 
sults in 0.5 < (M /Mn)~ 3 / 7 R fi 12 / 7 1 4Fl ~ 



: 1.6 



Wiri nga et al. 1988; lArnett and Bowers! 11977: 
Pandharipande et al. 19761) where Rk is in units 
of 10 6 cm and J45 is in units of 10 45 g cm 2 . The 
model of the energy spectrum was used to infer 
the bolometric fraction of the luminosity from 5 

- 50 keV. For OAO 1657-415 it was determined 
that 0.5 < a < 0.8. Although a cyclotron reso- 
nance scattering feature (cyclotron line) has never 
been detected for OAO 1657-415, the energy of 
the cyclotron lines in other X-ray binaries range 
between 10 - 60 keV. Broadening this range to 10 

- 100 keV implies a magnetic field strength range 
for this population to be 1.0 < B 12 < 10.0 G. The 
factor, (n(wo)A:) 1 / 2 , is assumed for this work, to 
be between 0.8 and 1.0. Using the fitted results in 
Figure [6] and Equation [Til the distance to OAO 
1657-415 can be constrained to between 3 and 1 6 
kpc which is consistent with iMason et al" (l2009h . 
4.4 < d < 12 kpc. The existence of a long lived 
(weeks to months) transient accretion disk is con- 
sistent with the flux-i> correlation and the distance 
to the star. 

A question of considerable importance to these 
type of systems is how a stellar wind can trans- 
fer nearly uniform (averaged over a few days) an- 
gular momentum to the neutron star for peri- 
ods of months. Even considering a stable wind 
emanating from t he companion which is unlikely 
(<Naik et al.ll2009l) . the disrupting presence of a 
neutron star and its accompanying X-rays irradi- 
ating the wind and the companion will add con- 
siderable complexity to the system. Aside from 
the short term periods of spin-up and spin-down, 




Fig. 7. — Models of the mass of a neu- 
tron star verses the combi nation of parameters , 



(M/Mp ) )~ 3 / 7 Rp, 12 ^ 7 hh~ 1 dWirinea et al.1 l"l 988: 



Arnett and Bowers! 119771 ; iPandharipande et al 
19761 ). The type of model is insensitive to this 



grouping of parameters. 



the long term decrease in spin period that has oc- 
curred since the discovery of the source suggests 
that the torque on the neutron star does not stem 
from stochastic processes. 

Although our temporal resolution for torque 
switching is a few days, Figure [6] is reminis- 
cent of an unstable prograde and retrograde 
transient (hours) accret ion disks, described by 
Taam and Frvxelll (|l989l) . forming around OAO 
1657-415 contributing to the torque responsible 
for v < 3 x 10~ 12 Hz s _1 . We have concluded 
that an accretion disk is present during the pe- 
riods in which v > 3 x 10~ 12 hz/s. These peri- 
ods last for weeks to months. It is possible that 
the disk is forming and disappearing on shorter 
timescales than this but our technique prevents us 
from determining torque-luminosity correlations 
on timescales shorter than the orbital period. 

Figure [S] shows the frequency, frequency deriva- 
tive, and the maximum of the orbital folded flux 
vs. time. During spin-up, OAO 1657-415 is 
more likely to exhibit torque-luminosity correla- 
tions which are not apparent while the source is 
spinning down. The correlated part of Figure [5] 
(v > 3 x 10~ 12 hz/s) spans the same flux range 
as the uncorrelated part, therefore the flux alone 
is not a good proxy for the torque applied to the 
neutron star (See the third panel in figure [8]) . If 
this is the case then the formation of the stable 
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26 f B6 
26,3* 




Fig. 8. — Top panel is a frequency history of OAO 
1657-415 showing, at the bottom of the panel, the 
intervals where the frequency derivative was calcu- 
lated. The intervals are travesed by a narrow light 
gray line at times where the frequency derivative 
and flux are not correlated. They are traversed by 
a wide dark gray line where the frequency deriva- 
tive and flux are correlated. The middle panel 
shows the calculated frequency derivative vs. time 
for the above intervals while the bottom panel 
shows the maximum of the Swift /BAT orbitally 
folded flux [10~ 10 erg cm~ 2 s _1 ] for the same in- 
tervals. 



prograde accretion disk is not triggered by high 
M but must have some other, unknown, trigger- 
ing mechanism. Additional three dimensional hy- 
drodynamic simulations are needed to understand 
these observations. 

Traditionally, the term wind-fed X-ray binaries 
has implied Bondi-Hoyle-Lyttleton accretion while 
disk-fed systems imply Roche-lobe overflow. Ob- 
servations of accreting X-ray binaries are begin- 
ning to depict a much more complicated picture 
with systems that include disk wind accretion and 
direct wind accretion systems. 
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Tabic 1: Orbital Ephcmcridcs for OAO 1657-415. 



Start MJD 


T T 


/ 2 a [MJD] 


a x sini a [lt- 


■sec] 


ecosw 




e sin a; 


x 2 


DOF 


48370. 


,7 


48390 


.6549±0. 


.0027 


106, 


,49±0 


,29 


-0.0053±0.0026 


0. 


1151±0. 


,0038 


9.8 


8 


48527. 


.5 


48547 


.3800±0. 


,0052 


104 


,01±0 


.61 


-0.0071±0.0027 


0. 


1120±0. 


,0046 


23.9 


8 


48558. 


.8 


48578 


,7293±0. 


.0050 


106 


,05±0 


.47 


-0.0011±0.0032 


0. 


1047±0 


.0050 


14.8 


8 


48715. 


.5 


48735 


,4386±0. 


.0033 


106 


,34±0 


,35 


-0.0032±0.0035 


0. 


1124±0. 


.0057 


7.6 


8 


49279. 


.7 


49299 


,5984±0. 


.0062 


105 


.09±0 


.60 


-0.0069±0.0037 


0. 


0999±0 


.0052 


23.8 


8 


49603. 


.6 


49623 


,4633±0. 


.0048 


106 


,42±0 


.52 


-0.0096±0.0034 


0. 


1101±0. 


.0052 


18.7 


8 


50240. 


.9 


50260 


,7701±0. 


,0026 


106 


,16±0 


.29 


-0.0053±0.0022 


0. 


1061±0. 


.0032 


12.3 


8 


50272. 


.2 


50292 


,1121±0. 


.0052 


106 


,82±0 


.56 


-0.0098±0.0036 


0. 


1140±0. 


.0056 


21.5 


8 


50303 


,5 


50323 


,4677±0. 


,0031 


106 


,90±0 


.34 


0.0061±0.0032 


0. 


1071±0. 


.0048 


1.5 


8 


50334. 


.9 


50354 


,8047±0. 


.0021 


105 


,40±0 


.23 


-0.0042±0.0021 


0. 


1047±0. 


.0032 


6.8 


8 


50564. 


.7 


50584 


,6562±0. 


,0045 


107 


,08±0 


.51 


-0.0035±0.0031 


0. 


1162±0. 


.0050 


16.4 


8 


54700. 


.0 


54721 


.7666±0, 


.0030 


106 


,30±0 


.30 


-0.0033±0.0027 


0. 


1045±0. 


.0048 


3.4 


8 


54732. 


,0 


54753 


.1092±0. 


,0034 


105 


,74±0 


.36 


-0.0086±0.0031 


0. 


1018±0. 


.0044 


12.2 


9 


55243. 


.0 


55254 


,5552±0. 


,0025 


105 


,84±0 


.30 


-0.0113±0.0024 


0. 


1015±0. 


.0036 


14.4 


9 


55286. 


.0 


55306 


,7930±0. 


,0025 


106 


,31±0 


,28 


-0.0056±0.0022 


0. 


1021±0. 


.0036 


14.8 


9 


55317. 


.0 


55338 


,1411±0. 


,0031 


107 


,00±0 


.44 


-0.0055±0.0024 


0. 


1114±0. 


.0038 


23.0 


9 


55505. 


.0 


55526 


.1813±0. 


,0024 


106 


,36±0 


.28 


0.0010±0.0025 


0. 


1098±0. 


.0040 


8.4 


9 


55765. 


.0 


55776 


.9135±0. 


,0038 


106 


,52±0 


.53 


-0.0064±0.0024 


0. 


1071±0 


,0036 


26.5 


6 



"Errors are statistical and inflated to reflect xi- 
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